Aging and Energy Landscapes: Application to Liquids and Glasses 



Stefano Mossa,^' ^' ^ Emilia La Nave,"'^ Francesco Sciortino,^ and Piero Tartaglia"'^ 

^ Dipartimento di Fisica and INFM Udr and Center for Statistical Mechanics and Complexity, 
Universita' di Roma "La Sapienza" Piazzale Aldo Mora 2, 1-00185, Roma, Italy 
'^Center for Polymer Studies and Department of Physics, 
Boston University, Boston, Massachusetts 02215 
(Dated: February 1, 2008) 

The equation of state for a liquid in equilibrium, written in the potential energy landscape for- 
malism, is generalized to describe out-of-equilibrium conditions. The hypothesis that during aging 
the system explores basins associated to equilibrium configurations is the key ingredient in the 
derivation. Theoretical predictions are successfully compared with data from molecular dynamics 
simulations of different aging processes, such as temperature and pressure jumps. 



Glasses are out-of-equilibrium (OOE) systems, char- 
acterized by slow dynamical evolution (aging) and his- 
tory dependent properties. Their thermodynamic and 
dynamic properties depend on the preparation method 
(cooling or compression schedules) as well as the time 
spent in the glass phase The extremely slow aging 
dynamics has often been considered an indication that 
a thermodynamic description of the glass state can be 
achieved by adding one or more history-dependent pa- 
rameters to the equation of state (EOS) [||, |[ |4[ |[|, |^, 

|, 10; Recent theoretical work — mostly based on 
mean field models of structural glasses, where analytic 
solutions of the OOE dynamics can be explicitly worked 
out — supports such possibility. 

In this Letter we derive an EOS for OOE conditions, 
based on a generalization of the potential energy land- 
scape (PEL) thermodynamic approach |4[ The 
EOS is derived under the hypothesis that during aging 
the system explores basins associated with equilibrium 
configurations, a condition which, as shown later, is sim- 
ple to implement in the PEL context. The proposed EOS 
for OOE conditions depends on one additional parame- 
ter, which for example can be chosen to be the average 
depth of the explored local minima of the PEL (the so- 
called inherent structures (IS) Q). The EOS allows 
for the first time a detailed comparison between predic- 
tions and "exact results" from out-of-equilibrium molec- 
ular dynamics simulations. We report such a comparison 
for a realistic model of the fragile liquid orthoterphenyl, 
one of the most studied glass forming liquids jl6j . 

The PEL formalism Q allows for a clear sepa- 
ration of vibrational and configurational contributions in 
thcrmodynamical quantities. Indeed, the equilibrium liq- 
uid free energy F{T, V) at temperature T and volume 
V is expressed as a sum of a vibrational contribution 
fvib^ representing the intra-basin free energy p^ , and a 
configurational contribution {—TSconf + &is)- Here the 
configurational entropy Sconf accounts for the number of 



explored basins of the PEL of average energy depth ejs- 
fyib, for a system of N atoms, depends on the curvature 
of the potential energy surface, i.e., on the 3A^ eigenvalues 
ujf of the Hessian matrix evaluated at the IS. Recently, 
calculating the V derivative of F{T, V), an analytic EOS 
has been derived, based only on the statistical properties 
of the landscape 

The separation in vibrational and configurational parts 
can be carried out for several properties. In our case, it 
is particularly useful to separate the instantaneous pres- 
sure as sum of two contributions: the pressure felt by the 
system in the inherent structure configuration {Pis), and 
an additional contribution related to the finite temper- 
ature, which we refer to as the vibrational contribution 



P{T, V, eis) = Pisieis, V) + F„,(r, V, as). 



(1) 
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In equilibrium ejs is a well defined function of T and V 
and, as well known, P depends only on T and V. Note 
that Pis, being the pressure experienced by the liquid 
in the IS, depends only on the basin's depth and on V. 
Instead, Pmb depends on the basin's depth as well as on 
T and V |§. 

For models for which the statistical properties of the 
PEL, i.e., the distribution of shapes and numbers of 
IS |l^ are known, Pjs and Pyib can be calculated 
theoretically. In other cases, Pis and Pvib can be calcu- 
lated numerically. In all cases, an landscape-based EOS 
based on Eq. (0) can be derived. 

The separation of the pressure in two contributions 
is crucial for the possibility of extending the equilibrium 
EOS to OOE conditions, since it provides a direct way to 
rebuild the pressure once T, V and e/g are known. In the 
PEL formalism, the hypothesis that the system explores 
basins associated with typical equilibrium configurations, 
i.e., configurations sampled in equilibrium conditions, is 
equivalent to assuming that the relations — calculated in 
equilibrium — linking Pis and Pyn, to e/5 and V are valid 
also in OOE conditions. If this is the case, the knowledge 
of eis, T and V is sufficient to calculate both Pis, Pvib 
and their sum P according to Eq. (|l]). Similarly, the 
values of P, T and eis are sufficient to predict V, since 
both Pis[eis,V) and Pyib{eis,T,V) can be estimated 
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FIG. 1: (a) Equilibrium relation between basin pressure Pis 
and basin depth e/s- Symbols are simulation results at the 
five volumes per molecule 0.337 nm'' (down triangles), 0.345 
nm'^ (circles), 0.353 nm^ (diamonds), 0.361 nm'' (up trian- 
gles), and 0.369 nm^ (squares) from Ref. Dashed lines are 
calculated using the parameterization discussed in Ref. 
Results at four intermediate volumes are also shown. The 
local minima are calculated using standard minimization al- 
gorithms [^. (b) Puii, as a function of T for the same five 
volumes. Each curve has been shifted by n x 20 MPa to avoid 
overlaps. Dashed lines show a suitable parameterization of the 
e/s, V and T dependence of P„,;6, as discussed in Ref. |17|. 



as a function of V . The predicted V is the value for 
which Pis{eis, V) + Pvib{eis, T, V) matches the external 
pressure. 

Next we apply the ideas discussed above to a system 
of iV = 343 molecules interacting via the Lewis and 
Wanstrom (LW) model for the fragile molecular liquid 
orthoterphenyl (OTP) The LW model is a rigid 

three-site model, with intermolecular site-site interac- 
tions modeled by the Lcnnard- Jones (LJ) potential [p2[ . 
We refer to Ref. |21 for all numerical details. An equilib- 
rium EOS, based on the assumption of a Gaussian land- 
scape, has been recently presented for this model [ pT[ , 
and successfully compared with the "exact" EOS calcu- 
lated from molecular dynamics simulations. 

Fig. |l| shows the relation between Pis and eis in equi- 
librium along different isochores and Pyib as a function of 
T for the LW model. Both quantities can be parameter- 
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FIG. 2: Paths in the P - T plane of the four OOE simula- 
tions discussed in the text. Arrows indicate: (i) a T-jump 
at constant V between liquid states; (ii) a P-jump at con- 
stant T between liquid states; (iii) a P-jump at constant T 
starting from a glass; (iv) a constant P heating of a glass. P 
and T-jumps have been simulated by changing the thermostat 
and barostat value at i = 0. The time constant of the ther- 
mostat and barostat has been fixed to 20 ps. The full lines 
indicate the equilibrium PiT) at five V values, from Ref. . 
The dashed lines connect the equilibrium state points and the 
glass states obtained by fast constant-1/ cooling. 



ized in a quite accurate way. Indeed, for this model, Pmh 
is weakly dependent on V and e/s and it is well repre- 
sented by a linear T-dependence. The parameterization 
of Pis and Pmh offers the possibility of estimating the 
system pressure once T, V and e/s are known. 

We study four different OOE protocols via computer 
simulation. The imposed external conditions are illus- 
trated in Fig. ||. We study i) a T-jump at constant V\ 
ii) a P-jump at constant T; iii) a P-jump at constant T 
in the glass phase and iv) an isobaric heating of a glass. 
For the cases i) and ii), conditions are chosen in such 
a way that within the simulation time window (50 ns) 
the system reaches the equilibrium state. In all stud- 
ied cases, averages over 50 independent simulations have 
been performed. 

Case (i). — Fig. ^ shows the comparison between the 
theoretical predictions and the numerical results for the 
constant- T-jump case. The e/5 values, calculated from 
the simulation (Fig. ||(a)), are used together with the val- 
ues of T and V as input to predict the evolution of P 
(Fig. p|(b)). Figs. p|(c) and (d) also show the comparison 

for Pis and W. W(e/s), defined as ln(a;i(e/s)), 
provides an indication of the similarity in basin shape 
(at the level of the harmonic approximation) between 
basins explored in equilibrium and in OOE conditions. 
The quality of the predictions indeed supports the valid- 
ity of the EOS and the hypothesis that the basins sam- 
pled during aging have the same relationship W{eis) as 
in equilibrium. 

Case (ii). — Fig. ^ shows the evolution of the system 
after a constant-T pressure jump. In this case, the e/s 
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FIG. 3: Dynamics after a T-jump, at volume per molecule 
V = 0.345 nm^ from 480 K to 340 K. (a) e/s; (b) P; (c) 
Pis', (d) W {tUo ~ Icm^^). Open symbols are simulation 
results, closed circles are theoretical predictions based on the 
OOE PEL-EOS and using as input the data shown in panel 
(a). Dashed lines indicate the known equilibrium values at the 
initial and final state. The error in the prediction at short and 
long times, when equilibrium conditions are met, provides a 
measure of the quality of the equilibrium EOS [jl7|, which 
agrees with the simulation data within ±5MPo. 



values calculated from the simulation are used together 
with the values of T and P as input to successfully pre- 
dict the evolution of V . For times shorter than 20 ps, 
the barostat time constant, the fast increasing external 
pressure forces the system to change the volume with a 
solid-like response. In this time window, the initial basins 
in configuration space are continuously deformed by the 
rapid_volume change, in agreement with the findings of 
Ref. 
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For time longer than 20 ps, when the exter- 
nal pressure has reached the equilibrium value, the sys- 
tem starts to explore basins different from the (deformed) 
original ones, and the time evolution becomes controlled 
by the aging kinetics. To support this interpretation we 
show in Fig. ^(e) the actual path in the ejs-Pis plane, 
where the change around 20 ps is clearly evidenced. 

Cases (in) and (iv). — The last two cases are respec- 
tively a constant T P-jump and a heating at constant P 
(with an heating rate of 1 K/ps), both starting from glass 
configurations. These initial configurations are generated 
by rapid constant-V^ quenches of equilibrium configura- 
tions, as shown schematically in Fig. ^. Fig. ^ shows 
the comparison between the theoretical predictions and 
the numerical calculations for V, using ejs as input. In 
the case of a P-jump (Fig. |^Top), again two different 
dynamical behaviors are observed: a fast dynamics pro- 
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FIG. 4: Dynamics after a P-jump at T = 320 K from 13.4 
MPa to 60.7 MPa . (a) e/s; (b) V; (c) Pis; (d) W. Open 
symbols are simulation results, closed circles are theoretical 
predictions based on the OOE PEL-EOS, using as input the 
data shown in panel (a). Horizontal dashed lines indicate the 
known equilibrium values at the initial and final state. The 
vertical dashed lines indicate the time at which the external 
pressure reaches the final value. Panel (e) shows the actual 
path of the aging process in the Pis-eis plane. 



cess describing the mechanical response of the glass to 
the external pressure change, followed by an extremely 
slow aging dynamics, during which basin changes are tak- 
ing place. In both cases and for all times, the out-of- 
equilibrium EOS is able to predict quantitatively the V 
changes. 

In summary, we have derived an extension of the PEL 
equation of state to model the thermodynamic properties 
of liquids under out-of-equilibrium conditions. Such an 
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FIG. 5: Top: Simulation of a P-jump from -270.1 MPa to 
— 115.4 MPa, starting from a glass configuration, obtained by 
quickly cooling to 70 K equilibrium configurations at volume 
per molecule V = 0.369 nm^ and T = 170 K. Bottom: Iso- 
baric heating, starting from the IS at V = 0.337 nm^ per 
molecule and T = 320 K. The heating rate is 1 K / ps. Cir- 
cles are the ejs values used as input to predict the evolution 
of V (lines). Squares are the V calculated directly from the 
simulations. 



extension — based on the hypothesis that an aging sys- 
tem evolves through regions of configuration space which 
are typically sampled in equihbrium — requires an addi- 
tional thermodynamic parameter which, for convenience, 
we have associated with the IS energy of the explored 
basins. 



The results reported in this Letter (Figs, ^-j^) show that 
the proposed generalization of the PEL equation of state 
is successful in predicting out-of-equilibrium thermody- 
namics, at least under the conditions and time scales 
probed by state-of-the-art numerical simulations. Larger 
T or P jumps and/or longer aging time could require 
a more detailed thermodynamic description, with more 
than one additional parameter. Despite such a possi- 
bility, the energy landscape approach developed in this 
Letter is a promising starting point for looking into more 
complex problems in the physics of OOE liquids and 
glasses. 
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